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Abstract 

The relative thermodynamic stability of surface reconstructions including vacancies, adatoms 
and additional layers on both CaO and Mn02 terminations is calculated to predict the surface 
phase diagram of (\/2 x -\/2)R45° CaMnOa (001) using ab initio thermodynamics. Stoichiometric 
and nonstoichiometric reconstructions are considered. A set of boundary conditions driven by 
the binary and ternary sub-phases from CaMnOa defines its bulk stability region, enclosing the 
stable surface reconstructions that can be in equilibrium with the bulk. Most of the surfaces take 
the magnetic ordering of the bulk ground state, G-type antiferromagnetic (AFM); however, some 
of the Mn02-terminated surfaces are more stable when the surface layer spins flip. At 573 K, 
the stoichiometric CaO- and Mn02-terminated surfaces are predicted to be thermodynamically 
stable, as well as a CaO-terminated surface reconstruction where half the Ca are replaced by Mn. 
The Mn02-terminated surface reconstructions dominate the phase diagrams at high temperatures, 
including phases with MnO and Mn02 adatoms per surface unit cell. 
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I. INTRODUCTION 


Perovskite CaMnOs has recently been the subject of increasing interest due to its strain- 
assisted multiferroicity 1, 2|, thermoelectric efficiency js, 4|, colossal magnetoresistance j^, 
and catalytic properties 6| . Its high catalytic activity for oxidation reactions and structural 
similarity to important biological catalysts has motivated studies on CaMnOa as a potential 
catalyst for water oxidation |7|. While atomic insight into the CaMnOa surface is imperative 
if its potential application as a catalyst is to be realized, very little has been investigated in 
this direction. Most of the previous studies have focused on magnetic properties |8l-lin| of 
the stoichiometric cubic CaMnOa surface. 

In this study, we determine the atomic and magnetic low-energy surface structures of or¬ 
thorhombic Pnma CaMnOa (001) as functions of environmental variables, developing surface 
phase diagrams by combining density functional theory (DFT) and thermodynamics 11 |. 
Various surface reconstructions for both CaO- and Mn02-terminated surfaces are explored 
within a (\/2 x \/2)i?45° surface periodicity. For structures with Mn atoms at the surface, 
we considered multiple magnetic spin conhgurations to ensure that the surface stability ac¬ 
counts for the potential magnetic reconstruction at the surface. The surface free energies 
are computed to generate the surface phase diagram within the bulk stability range. To 
illustrate the surface stability at a range of temperatures, the vibrational free energy is in¬ 
cluded in the calculation of the bulk stability region. Through our surface examination, we 
predict thermodynamically stable surface reconstructions at a range of conditions. 


II. METHODOLOGY 


A. Computational Methods 


DFT calculations were performed using the PBEsol 12|] generalized gradient approxi¬ 
mation (GGA), as implemented in the Quantum Espresso 13| computer code. Perdew 
and collaborators have shown that the PBEsol functional describes solids and surfaces ac¬ 
curately [^|. PBEsol-|- U was also used to test electronic properties, with Uf.s = 7.23 eV, in 
line with previous work [l^ and determined with the linear response scheme 15|. However, 


adding strong correlations to GaMnOs provides unsatisfactory description of the magnetic 
properties, as has been established by Luo et al. [l^. Spin-polarized electronic densities 
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were used in all calculations, treating the magnetic moments as collinear in all cases, with 
additional non-collinear (NC) calculations for some Mn-terminated surfaces and for the 
ternary s ub-p hases. All atoms are represented by norm-conserving, optimized , designed 


nonlocal 


17( 1 pseudopotentials with spin-orbit interaction (for NC magnetism). The pseu- 

treating the 3s, 3p, 3d, 4s, and 


dopotentials are generated with the opium package jl^ . 

Ap of Ca, the 3s, 3p, 3d, 4s, and 4p of Mn, and the 2s and 2p of O as valence states. An 
accurate description of the magnetic properties of all the systems studied is achieved when 
nonlinear core-valence interaction is included in the Mn pseudopotential by the partial core 


correction scheme 


19l. |20| . The CaMnOa magnetic ground state is collinear G-type anti¬ 


ferromagnetic (AFM) with an observed Neel temperature of about 130 K 2l|, estimating 


the magnetic interaction energy to be around 6.6 meV 


22j. Therefore, all calculations are 


run with a 70 Ry plane-wave energy cutoff to ensure accuracy for the small relative energies 
among the different magnetic conhgurations. The Brillouin zone is sampled using a 4 x 4 x 1 
Monkhorst-Pack 23| fc-point mesh for surface structures based on energy convergence for 
bulk CaMnOs. For binary and ternary sub-phases, a dense enough set of /c-points was used 
so that the total energy is converged. The ionic relaxation parameters are chosen so that 
the forces in the surface and bulk structures are lower than 2 meV/A. 


We use density functional perturbation theory 24, 


25( 1 to calculate phonon frequencies 


and vibrational displacement vectors at the F point for all binary and ternary sub-phases 
in order to include the vibrational free energy contribution into the dehnition of the free 
energy. Translational symmetry is ensured by enforcing the three translational acoustic sum 
rules requiring that the acoustic modes have zero frequency. This is accomplished through 
small corrections to the diagonal elements of the dynamical matrix. 


B. CaMnOs (001) Surface Structures 


The (001) surface of CaMnOs is the lowest energy, since it consists of alternating CaO 
and Mn02 planes, each of which is neutral in an ionic picture 26(]. We construct surface 
structures with different compositions and reconstructions for both CaO and Mn02 termi¬ 
nations, with symmetrical slabs of 7-9 layers and a vacuum gap between slabs of ~15 A. All 
atoms are fully relaxed. 

Over 80 surface terminations are made by varying the stoichiometry of Ca, Mn and O 
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TABLE I. Energies per formula unit and magnetic moments on Mn atoms for bulk CaMnOs in 
different magnetic phases. 



FM 

A-AFM 

C-AFM 

G-AFM 

Energy (meV/formula unit) 

38 

12 

4 

0 

Magnetic moment (ps/Mn) 

2.75 

2.68 

2.64 

2.61 


including additional Ca 3 ;Oy and Mn^-Oy layers in the (\/2 x \/2)i?45° surface symmetry. 
The positions of adatoms and additional layers are strategically selected; Ca and Mn atoms 
are placed on top of O atoms, and O species on top of Ca and Mn atoms. Each surface 
termination is identihed based on its stoichiometric termination (either CaO or Mn02) 
and whether its composition involves vacancies (-) or adatoms (+). As an illustration, 
Ca0+2.00 refers to a CaO-terminated surface with two oxygen adatoms per primitive unit 
cell of CaO. The number of adatoms and/or vacancies is tabulated per formula unit of 
either CaO or Mn02, even though there are two formula units in the surface super cell. 
Similarly, Mn02-0.5Mn-0.50 refers to a Mn02-terminated surface with one Mn vacancy and 
one oxygen vacancy per two Mn02 units. Illustrations of the thermodynamically stable 
surface terminations after relaxation are shown in Fig. [H These surface phases include (a,b) 
the stoichiometric surfaces, (c,d) their full surface coverage with adsorbed oxygen molecules, 
(e) the Mn02-terminated surface with one Mn vacancy and one oxygen vacancy, as well as (f) 
a CaO-terminated surface where a Mn adatom occupies a Ca vacancy. The last two surface 
phases (g,h) are stable at high temperatures but represent very exotic surface structures, 
as the “Mn adatoms and “O adatoms have different coordination than the stoichiometric 
atoms, potentially indicating active sites. 


The ferromagnetic (FM) and three AFM phases of bulk CaMnOs were analyzed. The 
A-AFM state has FM planes of alternating direction along the 2 ;-axis. The C-AFM arrange¬ 
ment consists of a Mn with nearest neighbors of opposite spin in the plane, and parallel spins 
between one plane and the next, while G-type AFM arranges the spins so that all nearest 
neighbors have opposite spin orientation. The magnetic ground state of bulk CaMnOs is 
found to be G-AFM Pnma insulator, as reported in previous theoretical and experimental 
studies 0,0 .a. Ou. computed e,e*o„.c band gap is 0.65 eV. s.gnificantiy unde.estinsat- 


ing the experimental value of 3.1 eV 


28| but in close agreement with previous computational 
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FIG. 1. Relaxed atomic structures of CaMnOs (001) thermodynamically stable surface re¬ 
constructions: Stoichiometric (a) CaO-terminated, (b) Mn02-terminated, (c) CaO-|-2.0O, (d) 
MnO2+2.0O, (e) Mn02-0.5Mn-0.50, (f) CaO-0.5Ca-|-0.5Mn, (g) Mn02-|-0.5Mn-|-1.00, and (h) 
MnO 2 -l-l.OMn-l-l.OO. Large light blue spheres: Ca; medium purple spheres: Mn; and small red 
spheres: O. All surface reconstructions are within a (\/2 x \/2)R45° surface periodicity, denoted 
with black dashed lines. Top and side views of the surface are shown for reconstructions on fig¬ 
ures (c-h). For (e), the dashed purple and red circles represent the Mn vacancy and O vacancy, 
respectively. For (g,h), the “O and “Mn atoms are O and Mn adatoms colored white and blue, 
respectively. 


work [271, l29|] The atom-projected magnetic moments of Mn atoms in each magnetic phase 
are shown in Table 2 as well as the relative energies with respect to the stable G-AFM phase, 
which are in very good agreement with those of a recent projector angmented wave calcn- 


lation nsing the PBE fnnctional 


29| . The snrface calcnlations are performed with G-AFM 
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ordering, magnetic spin flip at the surface and other possible magnetic spin conhgurations for 
surfaces with complex Mn 3 ;Oj/ reconstructions (See Supplementary Material). The magnetic 
spin flip conhguration is studied for all Mn02-based surfaces since stoichiometric Mn02 is 
found to be more energetically favorable when all the surface spins flip (making the top and 
second layers FM aligned) jsj. 


C. Thermodynamic Stability 


The free energies of CaMnOa (001) surfaces are calculated with the thermodynamic ap¬ 
proach described in this section. In this model, the CaMnOa bulk is considered to be in 
contact with an atmosphere composed of all its components in equilibrium. This implies 
that the system acts as a reservoir where removal and/or addition of any of its components 
at the surface leads to a range of thermodynamically stable reconstructions and composi¬ 
tions at different conditions. The oxygen chemical potential, poi is used to correlate the 
surface stability with environmental conditions such as temperature, T, and oxygen partial 
pressure, p( 02 ). 

The surface free energy, fl*, of an individual surface slab, i, is dehned as the excess amount 
of free energy needed to create the surface from its bulk form ll| . 


= A 

2A 


L7„ 


slab 




2A 


“ ^CahCa ~ -^Mrihlln ~ ^oho] ' 


( 1 ) 


The Gibbs free energy of the slab is denoted by The term represents the free 

j 

energy of the material in slab i as the sum of the chemical potentials from the bulk of each 
component, fij, multiplied by the number of its atoms in each slab, 30|]. The factor of 
i is included to account for top and bottom slabs. The surface free energy is reported as 
energy per unit area. 

The ^ term shown in Eq. (2) accounts for the off-stoichiometric atoms of component 
n with respect to component m in any surface slab i by relating the number of atoms of 
components and in the slab to their ratio in the bulk. 


T® = — 
2A 


jybulk \ 

]V® - iv® —— 

n yybulk j 

’'m / 
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( 2 ) 









In our case, the non-stoichiometric surface components were defined with respect to Mn 
atoms. Dehning with respect to Ca atoms does not change the results. Next, Eqs. (1) 
and (2) are merged to express the Gibbs surface free energy as 

^ [G'Lab “ ^Mn(hCa + hMn + 3/io)] 

~rMn,CahCa “ rMn,ohO- (3) 


Since the system acts as a thermodynamic reservoir, the chemical potentials of CaMnOs 
components are not independent; they are related to the chemical potential of CaMnOs 
crystal, /icaMnOs- The Gibbs free energy of bulk CaMnOs per formula unit, ^'caMnOs) then 
equal to pcaMnOs because the slabs must be in equilibrium with the bulk as shown below: 

/^Ca + /iMn + 3po = /^CaMnOs = fi'CaMnOg- (4) 


Inserting the relationship above into Eq. (3) simplihes the surface free energy to 


^ [l^Lab ^MnS'CaMnOs] 
“TMn.Cah-Ca “ TJ^^^ q/Tq- 


(5) 


The bulk stability region is dehned by a set of boundary conditions 8ll-l84l|. First, Ca 
and Mn are not allowed to form metallic precipitates by satisfying the following conditions: 


II <r 

hCa S fi'ca ) 


hMn < gUn- 


( 6 ) 

(7) 


Second, the bulk is stable, while binary metal oxides derived from its components do not 
precipitate: 


hCa + ho < gcM, (8) 

hMn + hO<^S> (9) 

3/iMn + 4/io < 5 'Mn 304 ) (^O) 

2hMn + 3hO ^ 5'Mn203) (^l) 

hMn + 2hO < 5'Mn02' (7^) 


Third, there should be no precipitation of other ternary sub-phases (CaajMuj^Oz) from the 
bulk CaMnOs. The experimental phase diagram of bulk CaMnOs illustrates the impor¬ 


tance of considering these ternary sub-phases, as they can co-exist with CaMnOs 




36|. 
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Consequently, the following conditions must be met: 


hCa + 2//Mn + 4//0 < 5 'Sn 204 ) 

(13) 

hCa + 7/rMn + 12/^0 ^ 5'CaMn70i2’ 

(14) 

2P'Ca + hMn + 4po < fi'Cai*Mn 04 > 

(15) 

3/^Ca + 2//Mn + 7po ^ fi'Ca 3 Mn 207 ) 

(16) 

4/iCa + 3/iMn + 10/ro < S'CaiMnsOio • 

(17) 


The DFT energy is taken as the Helmholtz free energy at 0 K with zero-point energy, so 
the Gibbs free energy and the DFT energy are related as follows ll| . 


Q ^ ^total ^ ^vib _ j^pV ^E + F 


conf 


-ivib 


(18) 


The energetic contributions provided by the pV and the conhgurational entropy 
terms to the t her mo dynamic stability are considered to be negligible, as has been proven in 


previous studies 
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34| 


The term which is equal to , is calculated from the vibrational modes of 

the system. The phonon density of states, cr(a;), is used to express as an integral over 
frequency, u\ 


F"*’- = / da;F"‘’-(r,a;)a(a;) 


(19) 


The analytical expression of the harmonic quantum vibrational free energy per mode as a 
function of temperature and frequency, F"’’(r,a;), can be written as 


kgTh-i (l - . 


( 20 ) 


Once the Gibbs free energy of each compound is calculated, the stability region of bulk 


CaMnOa as a function of T can be evaluated 


3l| . The vibrational 


for all bulk phases while it is neglected for all slab structures 


32 


free energy is included 


34| . The surface vibra¬ 


tional behavior is complex, and the vibrational modes could be affected by the electronic 
redistribution at the surface region; however, it has been shown for other systems that the 
phonon frequencies do not vary the order of stable surface phases [31 1- 

The chemical potential for each component of CaMnOs is dehned as its difference from 
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TABLE II. Experimental and computed formation energies (eV/formula unit) for binary and 
ternary sub-phases defining the bulk stability region of CaMnOs. The calculated formation energies 


for all Mn-based phases were performed with their corresponding magnetic ground-state structure, 


noted next to the compound. The experimental values are obtained from Ref. 





Experimental 

Theoretical 

CaO 

-6.59 

-6.71 

MnO (AFM-ID) 

-3.99 

-4.06 

Mn304 (EiM-IlF) 

-14.38 

-15.46 

MnaOa (EM) 

-9.93 

-11.36 

Mn02 (AEM) 

-5.39 

-6.60 

CaMn204 (AEM) 

— 

-18.23 

CaMn70i2 (NC-AFM) 

— 

-47.01 

Ca2Mn04 (G-AFM) 

— 

-20.56 

Ca3Mn207 (G-AFM) 

— 

-34.34 

Ca4Mn30io (G-AFM) 

— 

-48.12 

CaMn03 (G-AFM) 

-12.83 

-13.79 


“ EM layers in the [Ill] plane and successive antiparallel layers. 

^ Eerrimagnetic with Mn atoms (ttiitt) the corresponding order as Ref. 



the total energy of its reference state: 


AA^Ca = hca - + Fc, 

A/iMn RMii 


ivibjbulk \ 


Pbulk I p 


vibjbulk \ 
Mn ) 


Aro — Ro 2 (-^02 + ^02 


vib,gas 


( 21 ) 

( 22 ) 

(23) 


For the metals Ca and Mn, the reference energy is the DFT total energy and vibrational free 
energy of the elemental bulk crystal. For oxygen, the reference is the energy of an O atom 


an isolated O 2 molecule. The is corrected so that = 2Eo+ 33] since 

known shortco ming of DFT is the overestimation of binding energies [^, particularly for 


double bonds 


38 


39|. 


Combining the chemical potential expressions, Eq. (21-23), and the approximation to the 
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Gibbs free energy for the slabs and CaMnOs bulk leads to the following expression for the 
Gibbs surface free energy, 


G* = — 


1 

2A 


■^slab 


An / T^bulk _|_ pvibjbulk 
I -^CaMnOs -^CaMnOs 




.(24) 


The same approach including the vibrational free energy is applied to Eq. (4), whose ex¬ 
pression is rewritten introducing the formation energy, as follows. 


A/iCa + A/iMn + 3A/io — AEj 


bulk pvibjbulk 

fjCaMnOa '^CaMnOs' 


(25) 


This equation is then used to rearrange all the bulk stability boundary inequalities, Eq. 
(6-17), so that only pca and /Uq are kept as variables. For instance, Eq. (15) becomes 


A/XCa + A^o < AESmuO^ - A^fcaMnOa 

I 7;^vib,bulk pvib,bulk 

+ -^Ca2Mn04 “ “^CaMnOs* 


(26) 


The results from the DFT formation energies, AE}?'^*'^, for all the binary and ternary sub¬ 
phases are compared to experimental quantities in Table 2. The formation energies are 


slightly overestimated, which is a well-known shortcoming of DFT 




The surface phase diagram dependence on T and p(02) provides insightful physical in¬ 
terpretations of the experimental range of conditions for the surface stability. The Apo can 
be directly related to T and p(02) by the ideal gas approximation, since at equilibrium the 
po in GaMnOs bulk is equal to the chemical potential of oxygen gas in the environment. 
This leads to 


Ai^o(T,p) = Afio (T,/) + 


(27) 


where Aun( T,p^)=^AGf^( T,p^). Thermodynamic data from the NIST-JANAF thermo¬ 
chemical tables 42|] are used to determine the values of Afio{T,p^) by selecting the refer¬ 
ence state for ^AG^{T,p^) extrapolated to r=0 K and p(02)=l atm. Then, the values 
obtained are plotted for a temperature range of 0-3000 K at different p(02) to correlate the 
stable surface phases with controlled environmental conditions. 
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III. RESULTS AND DISCUSSION 


A. Magnetic Surface Reconstruction 


Depending on their snrface composition, some of the Mn02-based slabs nndergo a mag¬ 
netic spin flip reconstrnction at the snrface 8j that is energetically favorable, leading to 
lower total energy and greater snrface stability. This magnetic spin flip reconstrnction refers 
to a snrface with G-AFM bnlk order, bnt where the Mn spins in the snrface layer align with 
the spins of the Mn in the second layer, generating a bilayer with C-type AFM alignment, 
as shown in Fig. [2]J^a,c). This occurs when the surface breaks the degeneracy of Cg and t 2 g 
states, as the Mn cage is changed from octahedral to square pyramidal. This leads to the 
partial occupation of the Cg (1^2 bands at the surface, which favors the FM alignment of spins 
between the surface and subsurface Mn 02 layers. 


In bulk CaMnOs, the superexchange magnetic interactions between the hlled 1)2g shells 
and the formally empty eg shells determine the G-AFM ordering. However, this G-AFM 
ordering does not extend to the (001) Mn02-terminated surfaces, since double-exchange 
governs the magnetic interaction between Mn subsurface and surface. The Mn 3d orbital- 
projected density of states (DOS) of both G-AFM and spin flip surface and subsurface layers 
of the Mn02-terminated stoichiometric surface are shown in Fig. |2](b,d). The surface states 
of both G-AFM and spin-flip show the broken degeneracy of Cg —)■ dj. 2 _y 2 , d ^2 and hg -t 
{dxzydyz), as expected for square pyramidal. The G-AFM subsurface states have degenerate 
Cg and t 2 g shells; however, the spin-flip subsurface Cg degeneracy is broken, as the d ^2 is 
favorably aligned with the d^^ from its neighboring surface state. The alignment of the 
Mn spins on surface and subsurface layers is induced by a double-exchange process, which 
is enabled by the occupation of the d^^ orbital, lowering the energy of Mn02-terminated 
surfaces with the spin flip arrangement. Our DFT computed stoichiometric Mn02-surface 
yields a spin flip structure as its magnetic ground state, with an energy lowering relative to 
G-AFM of 15 meV per primitive surface cell. 


Based on the prevalence of the energetically favorable magnetic spin flip reconstruction on 
the stoichiometric Mn02-terminated surface, it is important to study the effect of magnetic 
reconstruction on surface stability. All the Mn02-based surface reconstructions were relaxed 
with both G-AFM and spin flip magnetic ordering. Different magnetic conhgurations were 
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FIG. 2. Spin-flip magnetic reconstruction at the stoichiometric Mn02-terminated surface. The 
atomic, electronic and magnetic structure of the stoichiometric Mn02-terminated surface adopting 
(a,b) the bulk magnetic configuration, G-AFM, and (c,d) an energetically more favorable magnetic 
reconstruction where the spins of the surface Mn atoms are flipped. The magnetic spins are 
represented by the yellow arrows. The orbital-projected densities of states (DOS) for surface and 
subsurface Mn d-orbitals of the (b) G-AFM, and (d) spin-flip magnetic structures show different 
energetic alignment of the 3d^2 orbital between surface and subsurface states. 


explored for CaO-terminated surfaces with Mn adatoms as well as Mn02 terminated surfaces 
with Mnj,Oy adatom reconstructions (See Supplementary Material). Half of the thermody¬ 
namically stable Mn02-based surfaces are magnetic spin flip conhgurations, as shown in Fig. 
E] (surface phases with red grid lines). For surface reconstructions with only one extra Mn 
adatom, the spin arrangement that makes a local FM alignment between the surface Mn 
adatom and the nearest Mn atom is more favorable. For example, the local FM alignment 
between the Mn surface adatom and the Mn on the subsurface layer in CaO-0.5Ca-|-0.5Mn, 


12 






















which is within the bulk stability region, is 12 meV per primitive surface cell more stable 
than a local AFM alignment (See Supplementary Material). 


B. Surface Phase Diagram 

The computed surface phase diagram of CaMnOs, generated after calculating the surface 
free energy of the atomic and magnetic lowest energy configuration of each surface phase, is 
shown in Fig. El This diagram predicts a total of fourteen thermodynamically stable surface 
phases in the given range of A/ica and A/iQ. The phase diagram can be divided into four 
different quadrants based on conditions of A/iq, A^uca and A/XMn- Even though A^Mn is not 
directly plotted, its value can always be deduced from Eq. (25). The upper left quadrant 
(I) shows oxygen rich and cation poor conditions; therefore, the favored phases are termina¬ 
tions with excess oxygen adatoms, denoted as MnO2-l-2.0O, Ca0-|-2.00 and CaO-|-1.0O. This 
quadrant also shows stability of the stoichiometric CaO- and Mn02-terminated surfaces as 
well as Mn02-based reconstructions, Mn02-0.5Mn-0.50 and MnO 2 -l-O. 5 Mn-l-l.OO. In quad¬ 
rant (II), both O and Ca are plentiful, leading to the absence of Mn02-terminated stable 
surfaces and the presence of Ca adatoms in stable CaO-based surfaces, CaO-|-1.0Ca-|-2.0O 
and CaO-|-1.0Ca-M.OO. The bottom right quadrant (III) involves oxygen-poor but relatively 
Ca-rich conditions. This region is characterized by stable surfaces with Mn02 termination, 
denoting that even though Ca is at high concentrations, Mn02-terminated surfaces are more 
favorable at oxygen poor conditions, clearly indicated by a surface phase where a Ca vacancy 
is hlled by a Mn adatom, CaO-0.5Ca-|-0.5Mn. In the left bottom quadrant (IV), Mn is at 
its richest condition since both O and Ca are at poor levels. This promotes the large stable 
region for both MnO 2 -l-l.OMn-l-l.OO and Mn02-l-2.0Mn-|-2.00. 

A boundary line from the top left to mid-bottom right corners of the diagram delimits 
the regions of CaO- (upper) and Mn02- (lower) terminated surfaces. Six different CaO- 
terminated surfaces are found, many of them covering a significant area of chemical potential 
space. On the other hand, Mn02-type phases generally show smaller stable regions with a 
higher variety of phases, which can be explained by the multivalency of Mn. For instance, 
the predicted path from Mn02-l-2.0Mn-|-2.00 to MnO2-l-2.0O consists of a complex series 
of transitions. These transitions occur by increasing the A/tq from -3.0 eV to -0.5 eV while 
keeping constant the A/ica at ~-7.0 eV. The changes in A/xq can be controlled experimen¬ 


ts 


tally by varying T or p{02)- For example, these surface phase transitions are predicted to 
occur at a constant p{02) of 10“^° atm, as the temperature is decreased from «il400 K to 
^400 K. Starting with the surface initially at ^1400 K in the Mn02+2.0Mn+2.00 phase, as 
the temperature is reduced to ~750 K, the desorption of oxygen atoms is thermodynamically 
favorable, leading to the surface reconstruction MnO2+1.0Mn+1.0O. Reducing the tempera¬ 
ture further to ?«600 K should cause the surface to reconstruct into MnO 2 -l-O. 5 Mn-l-l.OO and 
even co-exist with stoichiometric Mn02-terminated surface. Once the temperature reaches 
^400 K, the surface phase transition to MnO2-l-2.0O is thermodynamically favored. 

The slopes of the boundaries between phases directly relate to the changes in numbers of 
surface atoms of each species. The Ca0-|-2.00 surface phase changes to CaO-|-1.0Ca-|-2.0O 
at a vertical boundary of Apca ~-5.75 eV via adsorption of Ca cations. Similarly, horizontal 
lines at Apo ;^-0.8 eV and ~-0.9 eV denotes its transition to Ca0-|-1.00 and stoichiometric 
CaO, respectively, by the removal of oxygen adatoms. Between these extremes, the slope 
of the CaO to CaO-|-1.0Ca-|-2.0O transition, as indicated in Fig. [31 is -i, suggesting the 
addition of 2 Ca cations and 4 oxygen atoms per cell since the transition line will be crossed 
by moving toward Ca- and 0-rich conditions. Taking into account the relationship shown 
in Eq. (25), the slope of the boundaries between surface phases can then be dehned as 
[ 3 b-l) y where a, b, and c represent the net change in number of Ca, Mn, and O atoms, 
respectively, between phases. For example, the removal of 2 Mn atoms and the adsorption 
of 2 oxygen atoms per cell results in the transition of surface phase Mn02-l-2.0Mn-|-2.00 
to MnO 2 -l-l.OMn-l-l.OO. This is in agreement with the boundary for this phase transition 
having slope of -1. 


Stability Region 


The surface phase diagram shows a series of stable surfaces, but only a few of them 
are within the bulk stability region. The ternary sub-phases illustrated in Eqs. 13-17 are 
included in the dehnition of the stability region since an extensive experimental examination 
of the phase relations of Ca-Mn-0 bulk systems [35|, |36| found the co-existence lines for 
the CaMnOa phase with these ternary sub-phases, and correspondingly narrower stability 
boundaries for CaMnOa. The boundary of bulk phase stability is dehned by Ca 2 Mn 04 , 
CaaMn 207 , and Ca 4 MnaOio on the top and Mn 20 a, Mn 02 , CaMn 204 , and CaMn 70 i 2 on 
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FIG. 3. Surface phase diagram of CaMnOs with Mn metallic, binary, and ternary sub-phase bulk 
stability boundaries at r=0 K denoted as follows (l)Mn, (2)CaO, (3)MnO, ( 4 )Mn 304 , ( 5 )Mn 203 , 
(6)Mn02, ( 7 )CaMn 204 , ( 8 )CaMn 70 i 2 , ( 9 )Ca 2 Mn 04 , (10)Ca3Mn2O7, and (ll)Ca 4 Mn 30 io. The 
enclosed area shaded by black grid lines represents the bulk stability region. The phase diagram is 
divided into four regions, with dashed white lines to group phases with thermodynamic stability 
at similar condition ranges of A//o, A/rca; and A/iMn- The color-coded legend for each of the 
surface phases is on the left panel. The phases with red grid lines represent surfaces with spin-flip 
magnetic reconstructions. The right panel shows the T and ^(© 2 ) dependence of A//o- 


the bottom (lines 9-11 and 5-8, respectively, shown in Fig. [3]). 

A confirmation of the importance of including the ternary sub-phases is the ability of 
this model to describe the decomposition of bulk CaMn 03 . Experimentally, decomposition 


of CaMn 03 occurs at ~1700 K and 0.21 atm O 2 


35| . Using the DFT-generated stability 


region computed at 0 K, this decomposition temperature is predicted to be ^^2100 K at 0.21 
atm by considering the decomposition conditions at the low Apo end of the bulk stability 
region. However, if boundaries 7-11 are neglected, then the stability region is limited by 
Mn02, Mn 203 and Mn 304 on the bottom and bordered by CaO on the top (lines 6, 5, 4, and 
2, respectively, in Fig. [2]), artificially increasing the computed decomposition temperature 
to over 2500 K at 0.21 atm. 

These ternary sub-phases include the Ruddlesden-Popper phases CaO(CaMn 03 )„, which 
are perovskite structures with double CaO layer (001) planes with n layers of CaMn 03 in be¬ 


tween 


43l |. Naturally, the surface reconstruction with CaO double layer, CaO-|-1.0Ca-t-1.0O, 


is found on the side where the Ruddlesden-Popper phases CaO(CaMn 03 )„ limit the stabil- 
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ity of the CaMnOa bulk. The CaMn 204 phase is a marokite structure where the 
MnOe octahedra share edges and corners. This network of octahedra leads to channels 
along the a-axis which provide Ca atoms with 8-coordinated sites. It has been found that 


this marokite-type phase shows potential photocatalytic activity 


45| . This bulk phase bor¬ 


ders the surface phase with a Mn02-terminated layer that has an additional MnO layer, 
Mn02+1.0Mn-|-1.00. CaMn 70 i 2 is a complex distorted double perovskite, whose chemical 
formula is formally written as (CaMnj]“^)(Mnjj“^Mn’''^)Oi 2 . This perovskite-derived phase 


consists of j of the T-sites occupied by Ca anc 


are arranged via charge ordering below 440 K 


for some range of chemical potentials 31 
surface of CaMnOs. 




the rest by while the 5-site charges 


46l |. All these ternary sub-phases are stable 


dehning realistic stability boundaries for the 


Temperature Dependence 


The bulk stability region described above is based on DFT calculations at 0 K enclosing 
the CaO-|-2.0O, CaO-|-1.0O, CaO, CaO-0.5Ca-|-0.5Mn, MnO2+2.0O, and Mn02-0.5Mn-0.50 
phases. A more accurate description of the stability region requires the inclusion of the 
vibrational free energy. This is signihcant, since it determines how the stability region 
changes as a function of temperature. 

The bulk stability region is recalculated incorporating the bulk vibrational free energy 
at 573 K, 873 K and 1173 K, shown in Fig. 01 These temperatures are selected since 
they are within the range of experimental temperature conditions used in the synthesis 
and pulsed laser deposition of the CaMnOs crystal and him, as well as in the annealing 
process for surface characterization techniques. These temperatures shift the stability region 
such that only Mn02-based surfaces are favored. As temperature increases, the stability 
region shifts towards Ca-poor and 0-poor conditions, in agreement with previous studies 
on other surfaces {m, 33 |. As the stability region moves, it encloses other phases that 
are predicted to be stable at those temperatures. At 573 K, the surface phases within 
the bulk stability region are the same as in Fig. [3l with the addition of the stoichiometric 
Mn02-terminated surface. Increasing the temperature by 300 K leads to the stability of 
Mn02+1.0Mn-|-1.00 surface phase. At 1173 K, the phases within the bulk stability region 
include the Mn02+0.5Mn-|-1.00 phase. The clear trend is that the surface structures with 


16 

















0 



Ca0+2.00 

Ca0+1.00 

CaO-0.5Ca+0.5Mn 

CaO 

CaO+1.0Ca+2.0O 

CaO+1.0Ca+1.0O 

MnO2+2.0O 




MnO2-O.5Mn-O.5O 

Mn02 

Mn02+0.5Mn+1.00 

Mn02+1.0Mn+1.00 

Mn02+2.0Mn+2.00 

MnO2+0.5Mn 

MnO2+1.0Mn 



APca (eV) 


FIG. 4. Temperature dependence of the bulk stability region within CaMnOa surface phase dia¬ 
gram. Stability limits at different temperatures 573 K (yellow), 873 K (orange) and 1173 K (dark 
red) are plotted, indicating the thermodynamic stability of Mn02-terminated surface reconstruc¬ 
tions. 


Mn 02 terminations are more favorable at higher temperatnres. 

Changes in the shape of the stability region reflect how binary and ternary snb-phase 
bonndaries shift as the vibrational free energy at varions temperatnres is inclnded. Overall, 
the bnlk stability region tends to shrink as the higher temperatnre contribntions are inclnded. 
This is in good agreement with previons reports, since it is experimentally observed that 


CaMnOs crystal decomposes at temperatnres higher than «il700 K 
bnlk stability region vanishes aronnd the same temperatnre. 


35j and onr theoretical 


IV. SUMMARY AND CONCLUSIONS 

The energies of series of CaO- and Mn02-terminated CaMnOs (001) snrfaces with dif¬ 
ferent combinations of vacancies, adatoms and additional layers are reported nsing ab initio 
thermodynamics to theoretically predict the snrface phase diagram. The predicted snrface 
phase diagram has snrface strnctnres for both terminations, with stability ranges specihed 
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by A/io and A/ica- The A/zq is related to T and p{02), linking surface stability with 
experimental conditions. 


The chemical potential region for which CaMnOs is the only bulk phase is bounded 
by coexistence lines with the Ruddlesden-Popper phases Ca 0 (CaMn 03 )„ (where n=l, 2, 
and 3), as well as Mn 203 , Mn02, CaMn 204 , and CaMn 70 i 2 . This shows the relevance of 
incorporating the ternary subphases from the experimental Ca-Mn-0 compositional phase 
diagram 35|, ISq . The surface phases that are favored within this bulk stability region change 
once the vibrational free energy at different temperatures is included. The stability region 
shifts to lower A/xca for higher T, leading to the predominance of complex Mn02-based 
surfaces at T above 800 K. Thermodynamically stable phases include potentially reactive 
surfaces, Mn02+1.0Mn+1.00, since it shows two unique Mn sites based on coordination. 


The Mn02-based surfaces show intriguing atomic and magnetic structures which are com¬ 
prehensively explored by computing the magnetic ground state and analyzing its thermo¬ 
dynamic stability. Some of the Mn02-terminated surfaces are energetically more favorable 
when their surface and subsurface spins are aligned. This occurs when the hlling of the 
dz 2 orbital causes a more favorable double-exchange interaction between the subsurface and 
surface Mn atoms. The atomic composition of the surface affects this magnetic reconstruc¬ 
tion; for some Mn 3 ;Oy reconstructions, the surface reconstructs, rearranging the magnetic 
interactions in the system. However, this magnetic spin flip reconstruction prevails for key 
stable surfaces such as stoichiometric Mn02-terminated and MnO2+2.0O. 


This study addresses structures with {'/2x\/2)RA5° surface symmetry, leading to ordered 
surfaces due to periodic boundary conditions. Admittedly, real surfaces are more complex 
and disordered; however, by taking into account the chemical and thermal contributions to 
the surface stability and exploring a wide range of possible surface structures, the predicted 
phase diagram provides an acceptable foundation for the analysis of experimental data on 
this surface. Consequently, this predicted surface phase diagram could lead to experimental 
studies performing low-energy electron diffraction and scanning tunneling microscopy to 
validate the stable surface reconstructions and further explore the surface evolution. 
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